The binding site distance test score: a robust method for the assessment of predicted protein binding sites

We propose a novel method for scoring the accuracy of protein binding site predictions – the Binding-site Distance Test (BDT) score. Recently, the Matthews Correlation Coefficient (MCC) has been used to evaluate binding site predictions, both by developers of new methods and by the assessors for the community wide prediction experiment – CASP8. Whilst being a rigorous scoring method, the MCC does not take into account the actual 3D location of the predicted residues from the observed binding site. Thus, an incorrectly predicted site that is nevertheless close to the observed binding site will obtain an identical score to the same number of nonbinding residues predicted at random. The MCC is somewhat affected by the subjectivity of determining observed binding residues and the ambiguity of choosing distance cutoffs. By contrast the BDT method produces continuous scores ranging between 0 and 1, relating to the distance between the predicted and observed residues. Residues predicted close to the binding site will score higher than those more distant, providing a better reflection of the true accuracy of predictions. The CASP8 function predictions were evaluated using both the MCC and BDT methods and the scores were compared. The BDT was found to strongly correlate with the MCC scores whilst also being less susceptible to the subjectivity of defining binding residues. We therefore suggest that this new simple score is a potentially more robust method for future evaluations of protein-ligand binding site predictions.

Interactions of decay-accelerating factor (DAF) with haemagglutinating human enteroviruses: utilizing variation in primate DAF to map virus binding sites

A cellular receptor for the haemagglutinating enteroviruses (HEV), and the protein that mediates haemagglutination, is the membrane complement regulatory protein decay accelerating factor (DAF; CD55). Although primate DAF is highly conserved, significant differences exist to enable cell lines derived from primates to be utilized for the characterization of the DAF binding phenotype of human enteroviruses. Thus, several distinct DAF-binding phenotypes of a selection of HEVs (viz. coxsackievirus A21 and echoviruses 6, 7, 11-13, 29) were identified from binding and infection assays using a panel of primate cells derived from human, orang-utan, African Green monkey and baboon tissues. These studies complement our recent determination of the crystal structure of SCR(34) of human DAF [Williams, P., Chaudhry, Y., Goodfellow, I. G., Billington, J., Powell, R., Spiller, O. B., Evans, D. J. & Lea, S. (2003). J Biol Chem 278, 10691-10696] and have enabled us to better map the regions of DAF with which enteroviruses interact and, in certain cases, predict specific virus-receptor contacts.

In silico identification and characterization of protein-ligand binding sites

Protein–ligand binding site prediction methods aim to predict, from amino acid sequence, protein–ligand interactions, putative ligands, and ligand binding site residues using either sequence information, structural information, or a combination of both. In silico characterization of protein–ligand interactions has become extremely important to help determine a protein’s functionality, as in vivo-based functional elucidation is unable to keep pace with the current growth of sequence databases. Additionally, in vitro biochemical functional elucidation is time-consuming, costly, and may not be feasible for large-scale analysis, such as drug discovery. Thus, in silico prediction of protein–ligand interactions must be utilized to aid in functional elucidation. Here, we briefly discuss protein function prediction, prediction of protein–ligand interactions, the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated EvaluatiOn (CAMEO) competitions, along with their role in shaping the field. We also discuss, in detail, our cutting-edge web-server method, FunFOLD for the structurally informed prediction of protein–ligand interactions. Furthermore, we provide a step-by-step guide on using the FunFOLD web server and FunFOLD3 downloadable application, along with some real world examples, where the FunFOLD methods have been used to aid functional elucidation.

Evidence of a possible role for Lys-gamma 3-MSH in the regulation of adipocyte function

Lys-gamma 3-MSH is a melanocortin peptide derived from the C-terminal of the 16 kDa fragment of POMC. The physiological role of Lys-gamma 3-MSH is unclear, although it has previously been shown that, although not directly steroidogenic, it can act to potentiate the steroidogenic response of adrenal cortical cells to ACTH. This synergistic effect appears to be correlated with an ability to increase the activity of hormone sensitive lipase (HSL) and therefore the rate of cholesterol ester hydrolysis. Ligand binding studies have suggested that high-affinity binding sites for Lys-gamma 3-MSH exist in the adrenal gland and a number of other rat tissues that express HSL, including adipose, skeletal muscle and testes. To investigate the hypothesis that Lys-gamma 3-MSH may play a wider role in cholesterol and lipid metabolism, we tested the effect of Lys-gamma 3-MSH on lipolysis, an HSL-mediated process, in 3T3-L1 adipocytes. In comparison with other melanocortin peptides, Lys-gamma 3-MSH was found to be a potent stimulator of lipolysis. It was also able to phosphorylate HSL at key serine residues and stimulate the hyper-phosphorylation of perilipin A. The receptor through which the lipolytic actions of Lys-gamma 3-MSH are being mediated is not clear. Attempts to characterise this receptor suggest that either the pharmacology of the melanocortin receptor 5 in 3T3-L1 adipocytes is different from that described when expressed in heterologous systems or the possibility that a further, as yet uncharacterised, receptor exists.

High-affinity small molecule-phospholipid complex formation: Binding of siramesine to phosphatidic acid

Siramesine (SRM) is a sigma-2 receptor agonist which has been recently shown to inhibit growth of cancer cells. Fluorescence spectroscopy experiments revealed two distinct binding sites for this drug in phospholipid membranes. More specifically, acidic phospholipids retain siramesine on the bilayer surface due to a high-affinity interaction, reaching saturation at an apparent 1:1 drug-acidic phospholipid stoichiometry, where after the drug penetrates into the hydrocarbon core of the membrane. This behavior was confirmed using Langmuir films. Of the anionic phospholipids, the highest affinity, comparable to the affinities for the binding of small molecule ligands to proteins, was measured for phosphatidic acid (PA, mole fraction Of X-PA = 0.2 in phosphatidylcholine vesicles), yielding a molecular partition coefficient of 240 +/- 80 x 10(6). An MD simulation on the siramesine:PA interaction was in agreement with the above data. Taking into account the key role of PA as a signaling molecule promoting cell growth our results suggest a new paradigm for the development of anticancer drugs, viz. design of small molecules specifically scavenging phospholipids involved in the signaling cascades controlling cell behavior.

Mass spectroscopic characterization of the coronavirus infectious bronchitis virus nucleoprotein and elucidation of the role of phosphorylation in RNA binding by using surface plasmon resonance

Phosphorylation of the coronavirus nucleoprotein (N protein) has been predicted to play a role in RNA binding. To investigate this hypothesis, we examined the kinetics of RNA binding between nonphosphorylated and phosphorylated infectious bronchitis virus N protein with nonviral and viral RNA by surface plasmon resonance (Biacore). Mass spectroscopic analysis of N protein identified phosphorylation sites that were proximal to RNA binding domains. Kinetic analysis, by surface plasmon resonance, indicated that nonphospborylated N protein bound with the same affinity to viral RNA as phosphorylated N protein. However, phosphorylated N protein bound to viral RNA with a higher binding affinity than nonviral RNA, suggesting that phosphorylation of N protein determined the recognition of virus RNA. The data also indicated that a known N protein binding site (involved in transcriptional regulation) consisting of a conserved core sequence present near the 5' end of the genome (in the leader sequence) functioned by promoting high association rates of N protein binding. Further analysis of the leader sequence indicated that the core element was not the only binding site for N protein and that other regions functioned to promote high-affinity binding.

Binding of pentagalloyl glucose to two globular proteins occurs via multiple surface sites

The interaction between pentagalloyl glucose (PGG) and two globular proteins, bovine serum albumin (BSA) and ribulose-1,5-bisphosphate carboxylase oxygenase (rubisco), was investigated by isothermal titration calorimetry (ITC). ITC data fit to a binding model consisting of two sets of multiple binding sites, which reveal similarities in the mode of binding of PGG to BSA and rubisco. In both cases, the interaction is characterized by a high number of binding sites, which suggests that binding occurs by a surface adsorption mechanism that leads to coating of the protein surface, which promotes aggregation and precipitation of the PGG-protein complex. This model was confirmed by turbidimetry analysis of the PGG-BSA interaction. Analysis of tryptophan fluorescence quenching during the interaction of PGG with BSA suggests that binding of PGG leads to some conformational changes that are energetically closer to the unfolded state of the BSA structure, because small red shifts in the resulting emission spectra were observed.

Size and molecular flexibility affect the binding of ellagitannins to bovine serum albumin

Binding to bovine serum albumin of monomeric (vescalagin and pedunculagin) and dimeric ellagitannins (roburin A, oenothein B, and gemin A) was investigated by isothermal titration calorimetry and fluorescence spectroscopy, which indicated two types of binding sites. Stronger and more specific sites exhibited affinity constants, K1, of 104–106 M–1 and stoichiometries, n1, of 2–13 and dominated at low tannin concentrations. Weaker and less-specific binding sites had K2 constants of 103–105 M–1 and stoichiometries, n2, of 16–30 and dominated at higher tannin concentrations. Binding to stronger sites appeared to be dependent on tannin flexibility and the presence of free galloyl groups. Positive entropies for all but gemin A indicated that hydrophobic interactions dominated during complexation. This was supported by an exponential relationship between the affinity, K1, and the modeled hydrophobic accessible surface area and by a linear relationship between K1 and the Stern–Volmer quenching constant, KSV.

Hydrolyzable tannin structures influence relative globular and random coil protein binding strengths

Binding parameters for the interactions of pentagalloyl glucose (PGG) and four hydrolyzable tannins (representing gallotannins and ellagitannins) with gelatin and bovine serum albumin (BSA) have been determined from isothermal titration calorimetry data. Equilibrium binding constants determined for the interaction of PGG and isolated mixtures of tara gallotannins and of sumac gallotannins with gelatin and BSA were of the same order of magnitude for each tannin (in the range of 10(4)-10(5) M-1 for stronger binding sites when using a binding model consisting of two sets of multiple binding sites). In contrast, isolated mixtures of chestnut ellagitannins and of myrabolan ellagitannins exhibited 3-4 orders of magnitude greater equilibrium binding constants for the interaction with gelatin (similar to 2 x 10(6) M-1) than for that with BSA (similar to 8 x 10(2) M-1). Binding stoichiometries revealed that the stronger binding sites on gelatin outnumbered those on BSA by a ratio of at least similar to 2:1 for all of the hydrolyzable tannins studied. Overall, the data revealed that relative binding constants for the interactions with gelatin and BSA are dependent on the structural flexibility of the tannin molecule.

EfeO-cupredoxins: major new members of the cupredoxin superfamily with roles in bacterial iron transport

The EfeUOB system of Escherichia coli is a tripartite, low pH, ferrous iron transporter. It resembles the high-affinity iron transporter (Ftr1p-Fet3p) of yeast in that EfeU is homologous to Ftr1p, an integral-membrane iron-permease. However, EfeUOB lacks an equivalent of the Fet3p component—the multicopper oxidase with three cupredoxin-like domains. EfeO and EfeB are periplasmic but their precise roles are unclear. EfeO consists primarily of a C-terminal peptidase-M75 domain with a conserved ‘HxxE’ motif potentially involved in metal binding. The smaller N-terminal domain (EfeO-N) is predicted to be cupredoxin (Cup) like, suggesting a previously unrecognised similarity between EfeO and Fet3p. Our structural modelling of the E. coli EfeO Cup domain identifies two potential metal-binding sites. Site I is predicted to bind Cu2+ using three conserved residues (C41 and 103, and E66) and M101. Of these, only one (C103) is conserved in classical cupredoxins where it also acts as a Cu ligand. Site II most probably binds Fe3+ and consists of four well conserved surface Glu residues. Phylogenetic analysis indicates that the EfeO-Cup domains form a novel Cup family, designated the ‘EfeO-Cup’ family. Structural modelling of two other representative EfeO-Cup domains indicates that different subfamilies employ distinct ligand sets at their proposed metal-binding sites. The ~100 efeO homologues in the bacterial sequence databases are all associated with various iron-transport related genes indicating a common role for EfeO-Cup proteins in iron transport, supporting a new copper-iron connection in biology.

Thermodynamic analysis of ferrous ion binding to Escherichia coli ferritin EcFtnA

Iron oxidation in the bacterial ferritin EcFtnA from Escherichia coli shows marked differences from its homologue human H-chain ferritin (HuHF). While the amino acid residues that constitute the dinuclear center in these proteins are highly conserved, EcFtnA has a third iron-binding site (C site) in close proximity to the dinuclear center that is seemingly responsible for these differences. Here, we describe the first thermodynamic study of Fe2+ binding to EcFtnA and its variants to determine the location of the primary ferrous ion-binding sites on the protein and to better understand the role of the third C site in iron binding. Isothermal titration calorimetric analyses of the wild-type protein reveal the presence of two main classes of binding sites in the pH range of 6.5-7.5, ascribed to Fe2+ binding, first at the A and then the B sites. Site-directed mutagenesis of ligands in the A, B, or C sites affects the apparent Fe2+-binding stoichiometries at the unaltered sites. The data imply some degree of inter- and intrasubunit negative cooperative interaction between sites. Unlike HuHF where only the A site initially binds Fe2+, both A and B sites in EcFtnA bind Fe2+, implying a role for the C site in influencing the binding of Fe2+ at the B site of the di-iron center of EcFtnA. The ITC equations describing a binding model for three classes of independent binding sites are reported here for the first time.

Development of novel brush-type chiral stationary phases based on terpenoid selectors: HPLC evaluation and theoretical investigation of enantioselective binding interactions

The terpenoid chiral selectors dehydroabietic acid, 12,14-dinitrodehydroabietic acid and friedelin have been covalently linked to silica gel yielding three chiral stationary phases CSP 1, CSP 2 and CSP 3, respectively. The enantiodiscriminating capability of each one of these phases was evaluated by HPLC with four families of chiral aromatic compounds composed of alcohols, amines, phenylalanine and tryptophan amino acid derivatives and beta-lactams. The CSP 3 phase, containing a selector with a large friedelane backbone is particularly suitable for resolving free alcohols and their derivatives bearing fluorine substituents, while CSP 2 with a dehydroabietic architecture is the only phase that efficiently discriminates 1, 1'-binaphthol atropisomers. CSP 3 also gives efficient resolution of the free amines. All three phases resolve well the racemates of N-trifluoracetyl and N-3,5-dinitrobenzoyl phenylalanine amino acid ester derivatives. Good enantioseparation of beta-lactams and N-benzoyl tryptophan amino acid derivatives was achieved on CSP 1. In order to understand the structural factors that govern the chiral molecular recognition ability of these phases, molecular dynamics simulations were carried out in the gas phase with binary diastereomeric complexes formed by the selectors of CSP 1 and CSP 2 and several amino acid derivatives. Decomposition of molecular mechanics energies shows that van der Waals interactions dominate the formation of the diastereomeric transient complexes while the electrostatic binding interactions are primarily responsible for the enantioselective binding of the (R)- and (S)-analytes. Analysis of the hydrogen bonds shows that electrostatic interactions are mainly associated with the formation of N-(HO)-O-...=C enantio selective hydrogen bonds between the amide binding sites from the selectors and the carbonyl groups of the analytes. The role of mobile phase polarity, a mixture of n-hexane and propan-2-ol in different ratios, was also evaluated through molecular dynamics simulations in explicit solvent. (c) 2006 Elsevier Ltd. All rights reserved.

Isothermal titration calorimetry study of epicatechin binding to serum albumin

The interaction of epicatechin with bovine serum albumin (BSA) was studied by isothermal titration calorimetry. The binding constant (K) and associated thermodynamic binding parameters (n, Delta H) were determined for the interaction at three solution concentrations of BSA using a binding model assuming independent binding sites. These data show weak non-covalent binding of epicatechin to BSA. The interaction energetics varied with BSA concentration in the calorimeter cell, suggesting that the binding of epicatechin induced BSA aggregation. The free energy (Delta G) remained constant within a range of 2 kJ mol(-1) and negative entropy was observed, indicating an enthalpy driven exothermic interaction. It is concluded that the non-covalent epicatechin-BSA complex is formed by hydrogen bonding. (c) 2006 Elsevier B.V. All rights reserved.

SARS coronavirus unique domain: three-domain molecular architecture in solution and RNA binding

Nonstructural protein 3 of the severe acute respiratory syndrome (SARS) coronavirus includes a "SARS-unique domain" (SUD) consisting of three globular domains separated by short linker peptide segments. This work reports NMR structure determinations of the C-terminal domain (SUD-C) and a two-domain construct (SUD-MC) containing the middle domain (SUD-M) and the C-terminal domain, and NMR data on the conformational states of the N-terminal domain (SUD-N) and the SUD-NM two-domain construct. Both SUD-N and SUD-NM are monomeric and globular in solution; in SUD-NM, there is high mobility in the two-residue interdomain linking sequence, with no preferred relative orientation of the two domains. SUD-C adopts a frataxin like fold and has structural similarity to DNA-binding domains of DNA-modifying enzymes. The structures of both SUD-M (previously determined) and SUD-C (from the present study) are maintained in SUD-MC, where the two domains are flexibly linked. Gel-shift experiments showed that both SUD-C and SUD-MC bind to single-stranded RNA and recognize purine bases more strongly than pyrimidine bases, whereby SUD-MC binds to a more restricted set of purine-containing RNA sequences than SUD-M. NMR chemical shift perturbation experiments with observations of (15)N-labeled proteins further resulted in delineation of RNA binding sites (i.e., in SUD-M, a positively charged surface area with a pronounced cavity, and in SUD-C, several residues of an anti-parallel beta-sheet). Overall, the present data provide evidence for molecular mechanisms involving the concerted actions of SUD-M and SUD-C, which result in specific RNA binding that might be unique to the SUD and, thus, to the SARS coronavirus.

Titanocene anticancer complexes and their binding mode of action to human serum albumin: a computational study

Due to the pivotal role played by human serum albumin (HSA) in the transport and cytotoxicity of titanocene complexes, a docking study has been performed on a selected set of titanocene complexes to aid in the current understanding of the potential mode of action of these titanocenes upon binding HSA. Analysis of the docking results has revealed potential binding at the known drug binding sites in HSA and has provided some explanation for the specificity and subsequent cytotoxicity of these titanocenes. Additionally, a new alternative binding site for these titanocenes has been postulated.